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We determine the size of the critical region of the superfluid transition in the BCS-BEC crossover 
of a three-dimensional fermion gas, using a renormalization-group approach to a bosonic theory of 
pairing fluctuations. For the unitary Fermi gas, we find a sizable critical region [Tq,Tq], of order 
T c , around the transition temperature T c with a pronounced asymmetry: Tq — T c \/\Tq — T c \ ~ 8. 
The critical region is strongly suppressed on the BCS side of the crossover but remains important 
on the BEC side. 
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I. INTRODUCTION 

Superfluidity has been well understood for a long 
time in two different limiting cases. The first one 
is the Bardeen-Cooper-Schrieffer (BCS) superfluidity of 
fermions [1], the second one is the Bose-Einstein con¬ 
densation (BEC) of bosons [2]. In the BCS limit, a 
weak attractive interaction in a highly degenerate sys¬ 
tem of fermions induces a pairing instability. Pairs form 
and condensate at the same temperature T c that is or¬ 
ders of magnitude smaller than the Fermi energy Ep. In 
the BEC limit, superfluidity arises from condensation of 
bosonic atoms into a single quantum state. The inter¬ 
nal structure of the bosonic atoms (made up of an even 
number of fermionic constituents) is irrelevant at tem¬ 
peratures of the order of the BEC temperature. 

While most of the superfluids discovered in the 20th 
century are firmly in one or the other limit, recently 
discovered materials such as high-T c superconductors re¬ 
quire one to understand superfluidity beyond the stan¬ 
dard paradigms. From an experimental point of view, 
this has become possible in cold atomic gases where Fes- 
hbach resonances allow us to tune the attractive inter¬ 
action between atoms and span the entire BCS-BEC 
crossover [3]. Deep on the BCS side of the resonance, 
where the s-wave scattering length a between fermionic 
atoms in two different hyperfine states is negative and 
satisfies l/kp\a\ 1 (with n = k F /3n 2 the density of 
particles), superfluidity is well described by BCS theory. 
Deep on the BEC side, where a > 0 and 1/kpa 1, 
superfluidity arises from BEC of tightly bound fermion 
pairs of size a. The fermion gas close to unitarity 
(|a| —> oo) is a paradigmatic example of strongly corre¬ 
lated systems that exhibits several remarkable properties 
(for reviews, see Refs. [3-7]). In particular, it shows the 
highest ratio T c /Ep ~ 0.15 ever observed in a fermionic 
superfluid (with Ep = k 2 F /2m = (3n 2 n) 2 ^ 3 /2m) [8]. (We 
set H = kp = 1 throughout the paper.) 

An unsolved issue relates to the behavior of the unitary 
Fermi gas in the normal (nonsuperfluid) phase. In the 
BCS limit, above the transition temperature the Fermi 
gas is a Fermi liquid. On the BEC side, the normal phase 
is a Bose gas of fermion-fermion dimers which dissociate 


only at the pairing temperature T*; when T c < T <C T* 
the fermionic spectral function exhibits a pseudogap at 
low energies and Fermi-liquid behavior is suppressed. 
The importance of pairing fluctuations and “preformed” 
(metastable) pairs at unitarity as well as the fate of pseu¬ 
dogap [9, 10] and Fermi-liquid behavior is still unset¬ 
tled [ 1, 12] . On the experimental side, the situation 
is unclear. Radiofrequency spectroscopy has been inter¬ 
preted as evidence for a pseudogap state [13] and the 
small shear viscosity suggests the absence of quasipar¬ 
ticles above T c [14]. On the other hand, Fermi-liquid 
behavior was observed in thermodynamics [15, 16] and 
spin transport properties [17] in contradiction with the 
existence of a pseudogap. 

In this paper, we study pairing fluctuations in the 
BCS-BEC crossover and determine the size of the crit¬ 
ical region about the superfluid transition where ther¬ 
modynamic quantities show critical behavior with criti¬ 
cal exponents corresponding to the Wilson-Fisher fixed 
point of the three-dimensional classical 0(2) model. In 
weak coupling superconductors, due to their large T = 0 
coherence length, the critical region is too small to be 
observed experimentally; its size scales with ( T c /Ep ) 4 . 
For high-T c superconductors, it can reach 10~ 2 . It has 
previously been pointed out that the size of the criti¬ 
cal region in the unitary Fermi gas should be of order 
unity according to the Ginzburg criterion [ ]. However, 

a precise determination of the critical region requires a 
good description of critical pairing fluctuations and is 
therefore beyond the reach of standard theoretical ap¬ 
proaches such as Nozieres-Schmitt-Rink (NSR) theory 
or T-matrix approximation [4, 19, 20]. The superfluid 
transition temperature of the unitary Fermi gas has been 
determined accurately by a diagrammatic determinant 
Monte Carlo method [21] but the size of the critical re¬ 
gion has not been estimated. We have used the non- 
perturbative renormalization group (NPRG) to obtain 
an estimate of the critical region of the unitary Fermi 
gas [22], Our approach differs from previous NPRG stud¬ 
ies (for reviews, see Refs. [23, 24]) insofar as we integrate 
out the fermions using a complex Hubbard-Stratonovich 
field and deal with a bosonic theory of pairing fluctua¬ 
tions (for a perturbative RG approach along the same 
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lines, see [25]). 

We compute the Ginzburg temperatures Tq and Tq 
of the unitary Fermi gas and find a sizable critical region 
[Tq , Tq] around the transition temperature T c , of the or¬ 
der ofT c , with a pronounced asymmetry: \Tq— T c \/\Tq — 
T c | ~ 8. The critical region in the BCS-BEC crossover is 
estimated from a Ginzburg criterion with the additional 
constraint that the NPRG results should be reproduced 
in the unitary limit. We find that the critical region is 
strongly suppressed on the BCS side of the crossover but 
remains important on the BEG side. 


II. NPRG APPROACH 

We consider a system of neutral fermions with two pos¬ 
sible hyperfine states and dispersion e q = q 2 /2m, and 
Hamiltonian 

H = J d 3 r| $1 (-^ ' Vv - J^t} 

"" T4 _ . (1) 
where the chemical potential p controls the density of 
particles. We take g > 0 so that the interaction between 
fermions is attractive. The s-wave scattering length is 
then defined by 

^ = -i+f (2) 

4xa g J | q |< AF (2;r) 3 2e q 

where the ultraviolet momentum Ap on the fermion dis¬ 
persion is introduced to regularize the zero-range interac¬ 
tion potential gS( r — r'). We set h = kp = 1 throughout 
the paper. 

The partition function Z = Tre~^ ff can be written 
as a functional integral over anticommuting Grassmann 
fields with (Euclidean) action 

d 3 r'^2'ip*d r 'i/) tT +H[ip*,ip\ j. (3) 

Decoupling the interaction term in the particle-particle 
channel using a complex Hubbard-Stratonovich field ip 
and integrating out the fermion field i />, we obtain the 
action 



= - [ dr f d 3 r\ip{x)\ 2 - Trln(£ 0 1 + <£), (4) 
9 Jo J 


where j3 = 1/T and x = (r, r). 

G 0 (Qi i) = 


-l/„ . \ _ ( lUn - £q + d 0 


0 


lUJ n + £q~ fl 


(5) 


and 


ip(x, x') = S(x — x 1 ) 


0 <p(: r) \ 

<p(x)* 0 J 


( 6 ) 


are 2x2 matrices in the Nambu space defined by the 
two-component spinor 4/ = Here Go(q_- iu n ) 

denotes the Fourier transform of Qq(x,x') with u n = 
(2 n + 1)7 tT [n integer) a fermionic Matsubara frequency. 

The action S[ip*,<p\ contains all information about 
pairing fluctuations and criticality. A saddle- 
point computation of the partition function Z = 
f V[ip*, reproduces the BCS mean-field the¬ 

ory. Including Gaussian fluctuations about the saddle- 
point solution yields the NSR approach [19]. To take into 
account fluctuations beyond the perturbative approach, 
we use the NPRG [26-28]. We add to the action the 
infrared regulator term 

srv.d = E V (Pi * w m)Rfc (Pi 'iw m )</j(p, iu} m ) (7) 

p ,UJ m 

(oj m = 2rm:T is a bosonic Matsubara frequency) in¬ 
dexed by a momentum scale k such that fluctuations are 
smoothly taken into account as k is lowered from a mi¬ 
croscopic scale A down to zero. 

The central quantity in the NPRG approach is the 
scale-dependent effective action (or Gibbs free energy) 

T fc [A*,A]= - In Z k [J\J]+ dr J d 3 r{J*A + A*J) 

-^ eg [A*,A] (8) 

defined as a (slightly modified) Legendre transform of the 
(scale-dependent) free energy — In Z k [J*. J] that includes 
the subtraction of S™ s [A*, A]. Here J is a complex ex¬ 
ternal source that couples linearly to the pairing field ip 
and A (a;) = {<p(x)) is the order parameter. The cut¬ 
off function Ra(p) at the microscopic scale A is chosen 
very large so that it suppresses all fluctuations and the 
mean-field theory becomes exact: Ta[A*, A] = S[A*, A], 
For a generic value of k. R k (p) suppresses fluctuations 
with |p|, |w m |/c < k but leaves unaffected those with 
|p|, |w m |/c > k where c is the velocity of low-energy ex¬ 
citations. The effective action of the original model (4) 
is given by Tfc = o provided that R k -o vanishes. The vari¬ 
ation of the effective action with k is governed by Wet- 
terich’s equation [29], 

d k T k [A* , A] = * TV | ^ (Tf [A*, A]+74)^}, (9) 

where T).' denotes the second-order functional derivative 
of Tfc. Thus the NPRG approach aims at finding r^ = o 
using Eq. (9) starting from the initial condition Ta- 

At momentum scale A, all physical quantities of in¬ 
terest can be obtained from Ta[A*,A] = 5[A*,A]. In 
particular, the state of the system (normal or super- 
fluid) is determined by the constant (i.e., uniform and 
time-independent) minimum Aa of Ta- At unitarity, Aa 
becomes nonzero below the mean-field transition tem¬ 
perature T^ F ~ 0 AEp. In the following, we want to 
study the vicinity of the (true) transition temperature 
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T c ~ 0.15 Ep which is much smaller than T^ F [21], 
This implies that Aa is a large energy scale of the or¬ 
der of Aa (T = 0) ~ 0.69 Ep [30]. Restricting ourselves 
to fluctuations satisfying |p|,|w|/ca < Aa (with ca the 
Bogoliubov-Anderson mode velocity at mean-field level), 
we can expand the effective action Ta to second order in 
derivatives. Assuming a similar derivative expansion for 
all values of k smaller than A, this leads us to the ansatz 

r fc [A*, A] = f^dT I d 3 r{A*(Z c , k d T -V A , k d 2 

-Z A , k V 2 )A + U k (p)} (10) 

for the scale-dependent effective action, where the ef¬ 
fective potential U k is a function of the U(l) invariant 
p = |A| 2 . We further expand U k to quartic order in the 
field about the equilibrium state p = po,fc defined by the 
minimum of U k , 

U k (p) = U k (p 0tk ) + S k (p — po.fc) + — Po, k ) 2 , (11) 


where S k = U k (po, k ) vanishes if po. k > 0. 

Since the derivative expansion (10) is appropriate 
only for low-energy fluctuations satisfying |p|, |w m |/cA < 
Aa/ca, the initial momentum cutoff A must be of the 
order of Aa/ca- We thus assume that high-energy fluc¬ 
tuations do not significantly alter the transition tempera¬ 
ture and the size of the critical region. The precise value 
of the cutoff A is fixed so as to reproduce the known 
transition temperature T c — 0.15 Ep of the unitary Fermi 
gas [21], and internal consistency of our approach requires 
A ~ Aa/ca- 

The initial values 5a and Aa are obtained from 


<5a 


dU A (p) 

dp 


P=PO,A 


Aa 


d 2 U A (p) 

dp 2 


P—P 0,A 


( 12 ) 


whereas 


Zc ,a = lim w^-r^ 2) 12 (p), 

P ~^0 OUJm 

y ^ A = iToa^ r S 2(p) ’ (13) 

Za ' K = lim J^2 T a&(p)- 


U a is the mean-field effective potential and 



,5 2 r A [A*,A] 

SA i (-p)6A j (p) 


A(x)=A a 


— + ^—^- f to\Q A (q)T l Q + p)t j ] 

9 * J q 

(14) 


the mean-field two-point vertex in a uniform and time- 
independent field A (see appendix). We use the no¬ 
tation p = (p ,iu> m ) and q = (q, ioj n ) with w m (u> n ) 


a bosonic (fermionic) Matsubara frequency. S A x (q) = 
&o 1 (q) + Aat 1 is the mean-field fermion propagator and 
(t 1 ,t 2 ,t 3 ) stands for the Pauli matrices. Ai and A 2 re¬ 
fer to the real and imaginary parts of the pairing field 
A = (Ai +*A2)/V2. With no loss of generality, we have 
assumed Aa = p\j / 1° b e rea l- 

The bosonic effective action (10) has been studied in 
the context of superfluidity in Bose gases and we refer to 
previous publications for details about the solution of the 
flow equation (9) [31-34]. We have used an exponential 
cutoff function, 

R k (p) = Z A , k Yr(Y), r(Y) = (e y - l)" 1 , 

p 2 ^ (15) 

k 2 k 2 c\ ’ 

acting both on momenta and frequencies [33]. 


III. CRITICAL REGION OF THE UNITARY 
FERMI GAS 

To study the unitary limit, we compute both the ef- 

( 2 ) 

fective potential U A and the two-point vertex T ^' in 
the limit Ap — > 00 and g —> 0 in such a way that 
1/a = 0 (see appendix). We then fix the chemical po¬ 
tential p to an arbitrary value and determine the su¬ 
perfluid transition temperature T c from the vanishing 
of po.fc for k —> 0. Finally we fine tune the value of 
the ultraviolet cutoff A so as to obtain T c = 0.15Ep 
where n = — dU k= o(po, k =o)/dp is the fermion density and 
Ep = kp/2m = (37r 2 n) 2 / 3 /2m. We find A ~ 0.97fcp 
to be of the order of Aa (T = 0)/ca = 0.77kp [30]. 
When we compute physical quantities (e.g., the corre¬ 
lation length) away from the transition temperature, we 
adjust the value of the chemical potential so as to main¬ 
tain the fermion density n constant. 

Figure 1 shows the superfluid correlation length 

<“> 

in the normal phase. Near T c , £ ~ (T — T c )~ u diverges 
with the critical exponent v — 0.62, which is the known 
result for the three-dimensional classical 0(2) model with 
the approximations used here to solve the RG equa¬ 
tion (9). Defining the Ginzburg temperature Tq by the 
criterion 


In (£k F ) - Infefcf) 
In (£, c k F ) 


(17) 


where £ c ~ t~ u corresponds to the asymptotic behavior of 
the correlation length for t —> 0, we find that the critical 
region [T c . X/t] extends up to the Ginzburg temperature 
Tq ~ 1.75T C ~ 0.26 Ep. The Ginzburg temperature is 
of course sensitive to the precise criterion that we use. 
For example, one finds Tq ~ 1.66T C (and Tq ~ 0.94T C ) 
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FIG. 1. (Color online) Superfluid correlation length £ vs t = 
( T—T c ) /T c in the normal phase (T > T c ). The (green) dashed 
line shows the critical behavior £ ~ t~ v when t —► 0 + . (Inset) 
Ratio between the correlation length £ and the fermion de 
Broglie thermal wavelength £ t h = vf/kT. 



FIG. 2. (Color online) Order parameter y^po vs t = (T — 
T c )/T c in the superfluid phase (T < T c ). The (green) dashed 
line shows the critical behavior po ~ (—t) ,s when t —> 0 _ . 


with 0.01 instead of 0.02 in (17), and Tq ~ 1.97T C (and 
Tq ~ 0.79T C ) with 0.05. 

Figure 1 (inset) also shows the ratio between the corre¬ 
lation length £ and the fermion de Broglie thermal length 
£th = vf/ttT ( vf = kp/m). £/£th is of order one for 
T > Tq but rapidly increases for temperatures below 

T + 

1 G' 

We show in Fig. 2 the order parameter y/po = y/po,/c=o 
in the superfluid phase. Near T c , y/po ~ (T c — T)P van¬ 
ishes with the exponent /3 ~ 0.33. The critical behavior 
extends from the transition temperature T c = 0.157?^ 
down to the Ginzburg temperature Tq ~ 0.90T C ~ 
OAAEp- The latter is obtained from the criterion (17) 
with t^kp replaced by y/ ~po/Ep and y/po,c ~ (— t)&. 

We conclude that the critical region [ Tq,Tq} is sizable 
and strongly asymmetric: (Tq — T C )/(T C — Tq) ~ 8. 
This asymmetry is partially explained by the fact that we 
vary the temperature at fixed density while the chemical 
potential p strongly varies with T in the normal phase. 
At fixed chemical potential we find a smaller asymmetry: 


(T+ - T C )/(T C - Tq) cs 3. 

To clarify the origin of this asymmetry we have studied 
the classical 0(2) model (<p 4 theory for a two-component 
field) using the NPRG. We have found that the asymme¬ 
try is small when the interaction is weak but can be large 
when the interaction is strong. This shows that the criti¬ 
cal regime is generically asymmetric when the interaction 
is of order one, which is the case in the unitary Fermi gas. 
However, in the 0(2) model the critical regime is larger in 
the low-temperature phase, i.e., Tq > Tq , while the re¬ 
verse is true in the unitary Fermi gas. A crucial difference 
between the 0(2) model and our theory of the unitary 
Fermi gas is that in the latter case, the initial condi¬ 
tion of the RG approach is given by the mean-field BCS 
theory and is temperature dependent. By contrast, the 
only temperature dependence in the 0(2) model comes 
from the coefficient of the quadratic (<p 2 ) term. This, we 
believe, might explain the difference between the 0(2) 
model and the unitary Fermi gas but we have not been 
able to find a simple physical explanation. 

IV. CRITICAL REGION IN THE BCS-BEC 
CROSSOVER 

A priori it should be possible to follow a similar proce¬ 
dure to compute the size of the critical regime in the 
whole BCS-BEC crossover. While this is indeed the 
case on the BCS side (see below), in the BEC limit we 
encounter the following difficulty. When 1/kpa 1, 
superfluidity is due to the BEC of pointlike composite 
bosons. The initial effective action then corresponds to a 
bosonic action with a local interaction. To leading order 
in n 1 / 3 a, the transition temperature is given by the BEC 
temperature of bosons with mass ms = 2 m and density 
rip = n/2 . As is well known in the theory of the dilute 
Bose gas, the BEC temperature is not correctly com¬ 
puted if one imposes a cutoff on frequencies. Technically 
this is due to the propagator of the bosonic field (i.e., 
A) vanishing as 1 /iu> m in the large frequency limit. The 
calculation of the particle density then requires one to 
sum over all Matsubara frequencies (with an appropriate 
convergence factor). We expect that a cutoff acting only 
on momenta would allow us to explore the BEC limit but 
one would then have to go beyond the derivative expan¬ 
sion (for the frequency dependence), a very difficult task 
in practice. 

To estimate the critical regime in the BCS-BEC 
crossover and in particular in the BEC limit, it is pos¬ 
sible to use the Ginzburg criterion, with the additional 
constraint that the NPRG results should be reproduced 
in the unitary limit. Let us consider the classical limit 

F cl [A*,A]=/3 J d 3 r{Z A \VA\ 2 + U(p)} (18) 

of the effective action Tfc = A obtained by considering 
only fluctuation modes with vanishing Matsubara fre¬ 
quencies [35]. We can then define two characteristic 
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FIG. 3. (Color online) Critical region of the superfluid transi¬ 
tion in the BCS-BEC crossover. The transition temperature 
T c is obtained from a Luttinger-Ward approach [11]. The 
Ginzburg temperatures Tq and Tq (dashed and dot-dashed 
lines) are deduced from (19). The black crosses on the BCS 
side show the results obtained directly from the NPRG ap¬ 
proach (as in the unitary limit). 


ing the procedure used in the unitary limit (the cutoff A 
is fixed in order to reproduce the transition temperature 
computed in Ref. [41]). There is a good agreement with 
the results obtained from the Ginzburg criterion. By 
contrast, on the BEC side 1/fc^a > 0, the NPRG ap¬ 
proach yields a much narrower critical region than what 
is predicted by the Ginzburg criterion. We ascribe this 
disagreement to the failure of our NPRG approach (with 
a frequency cutoff) to determine correctly the particle 
density when there are composite bosons (bound pairs 
of fermions); see the discussion at the beginning of this 
section. 

The overall picture of the critical region is in qualita¬ 
tive agreement with the results of Ref. [18]. Evidence for 
a large critical region close to unitarity was also found 
in Ref. [42], where a diagrammatic analysis showed that 
diagrams that are subleading for T>T C exceed leading 
ones in a wide region, of order T c , above the transition 
temperature. 


lengths. The first one is the T = 0 healing length 
£h = ( Za/Xpo) 1//2 ; the second one is the Ginzburg length 
= Z\/\T C [36]. The size of the critical region is de¬ 
termined by the Ginzburg criterion 


‘'G c2 ’ 


T r 




(19) 


where a± is a constant. Here we distinguish between 
the critical region in the normal (t^) and superfluid ( t q) 
phases. Gi. and £g are computed vs 1 /kpa using mean- 
field theory (see appendix) and we fix the constants a± 
by requiring that for 1/a = 0 we reproduce the NPRG 
results obtained for the unitary gas in Sec. Ill: a+ ~ 
4.0 x 10~ 3 and a_ ~ 5.2 x 10 -4 . The asymmetry of the 
critical region is then independent of 1 /fcpa and equal 
to a+/a_. We therefore do not expect the Ginzburg 
criterion (19) to be quantitatively valid beyond the BCS- 
BEC crossover region l/kp\a\ < 1. 

Using the results of the appendix, in the BCS limit we 
find that t^, ~ ( T c /Ep) 4 is exponentially small, since 
T c ~ Epe - n / kF \ a \, so that the critical region is too 
small to be observed experimentally. In the BEC limit, 
t/i ~ kpa(T c /Ep) 2 ~ kpa [18, 37] is also small although 
not as strongly suppressed as in the BCS limit. For 
1 Ikpa 1, using T c ~ 0.218.EF (and the results of the 
appendix) [38], we obtain Tq ~ 7n 4 / 3 a <C 1 in agreement 
with known results of the dilute Bose gas [39, 40], By 
contrast, t q is of order unity at unitarity. Note that the 
maximum of Tq is on the BEC side of the resonance 
(1 /kpa — 1), not far away from the maximum of T c . 

The critical region in the BCS-BEC crossover is shown 
in Fig. 3, where we have used the transition tempera¬ 
ture T c computed in Ref. [41] within a Luttinger-Ward 
approach. On the BCS side 1 /kpa < 0, we also show the 
critical region obtained from the NPRG approach follow- 


V. CONCLUSION 

In summary, we have determined the width of the 
critical region in a Fermi gas through the BCS-BEC 
crossover. To circumvent the breakdown of perturba¬ 
tion theory in the critical region, we have used a NPRG 
approach. Close to unitarity, we find that the critical 
region is very wide (of the order of the transition tem¬ 
perature T c ) and asymmetric, being larger in the normal 
phase than in the superfluid phase. The critical region is 
also more pronounced on the BEC side than on the BCS 
side where it becomes unobservable. An interesting and 
open issue relates to the interplay between critical fluc¬ 
tuations, Fermi-liquid behavior, and the (much debated) 
possible existence of a pseudogap [43]. 

In Ref. [18], it was pointed out that the critical region 
is strongly suppressed in trapped gases due to their inho¬ 
mogeneity, although critical behavior could be observed 
in the compressibility by measuring the density profile. 
Nonetheless, the possibility to load gases in the three- 
dimensional quasiuniform potential of an optical box trap 
opens up new prospects to observe criticality in trapped 
gases [44]. 
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APPENDIX: EFFECTIVE POTENTIAL AND 
TWO-POINT VERTEX IN THE MEAN-FIELD 
APPROXIMATION 


At the initial stage of the RG procedure, the effec¬ 
tive action is given by mean-field theory: Ta[A*,A] = 
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S[A*, A], where S is given by (4). An elementary calcu- We use the notation 
lation then gives 


U A (p,p) = - + [ U q - 2Tln[2cosh(/3£ q /2)]} (20) 

9 J q 

d 

r A,ii(p; p) = - + \ [ { 2F (q) F (p + q) 

9 ^ J q 

- G(q)[G(p -q) + G(-p - 9)]}, 

r A,22(p; p) = --\t {2 F (Q)P(p + q) (21) 

9 Z J q 

+ G(q)[G(p -q) + G(—p - 9)]}, 

r A,i2(p; p)= G (q)i G (p -q)~ g (-p - 9)], 




d 3 q 
(2t r) 3 ' 


(23) 


The fermionic Matsubara sums in (21) can be done an¬ 
alytically. The interaction constant g can be eliminated 
using (2). We then take the limit Ap —> oo in the mo¬ 
mentum integrals with a fixed. 


in a constant field A = yfp (i.e., Ai = y/2p and A 2 = 0), 
with 


G(q) 


+ £ q . . _ y/p 

ul + EX W 0,2+i^’ 



( 22 ) 


Derivative expansion 

Equations (20) and (21) yield 


Aa = 


^3 tanh(/3E q /2) - ^ cosh 2 (/3£ q /2) 


ZA ’ A - i I 


2q 2 

dm 


^3 tanh(/3£ q /2) - ^ cosh 2 {/3Ej 2) 


8 m 2 d 


3/3 


q 2 e 2 is tanh(/3E q /2) - ^ cosh- 2 (^E q /2) - 




2E4 


/3 2 tanh(/3i? q /2) 
2E|cosh 2 (/3E q /2)J 


Zc, A = )/?, 


^3 tanh(/3E q /2) - ^ cosh 2 (/3£ q /2) 


VA,a = 


^cosh-^,/ 2 )- 


(24) 


where d = 3 and the right-hand side is evaluated at the 
minimum of the effective potential defined by [/)((/9o,a) = 
0. 

In the T = 0 BCS limit, we obtain 


U\{p) ^ -©(/*) j^(2w) 3 /2 M 5/2 , 


and 


VPo.a = 8e 2 e F e A a = 


Za, a = 


Z(ja — 0, Va,a = 


mkF 
47T 2 p 0 ,A ’ 
mfcp 


(25) 


(26) 


8mp o y 87T 2 p 0 ,A ’ 

where n = — dU\/dp = k 3 F /37 t 2 , fc F = y/2mtF 1 and e F = 

p. 


In the T = 0 BEC limit, we find 


mp p(2m) 3 / 2 |/Lt| 1 / 2 p 2 (2m) 3 / 2 |/r| 3 / 2 

t/A( " ) = - to + -8i- + -256i- 


and 


47m 


P 0,A — —2 


m*a ma* 


j5p, Aa = 


ma m 2 a 

^’ A= m’ Zc ’ A= ^’ 


167T 

Va,a = 


to 3 a 3 


(27) 


( 28 ) 


32tt 


where p, = — l/2ma 2 + 5/x (0 < (5/7 <C l/2ma 2 ). 
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